CURRENT AND DENSITY FLUCTUATIONS FOR INTERACTING 
PARTICLE SYSTEMS WITH ANOMALOUS DIFFUSIVE BEHAVIOR 



M. JARA 

Abstract. We prove density and current fluctuations for two examples of 
symmetric, interacting particle systems with anomalous diffusive behavior: the 
zero-range process with long jumps and the zero-range process with degenerated 
bond disorder. As an application, we obtain subdiffusive behavior of a tagged 
particle in a simple exclusion process with variable diffusion coefficient. 



1. Introduction 

Since the work of Harris [15], it is known that the motion of a tagged or dis- 
tinguished particle in symmetric, diffusive, one-dymensional systems of particles 
that preserve the relative order of particles is subdiffusive, in the sense that the 
mean square displacement E[xQ{t)'^] grows with t as t^^^, where denotes the 
position of a tagged particle, initially at the origin. This is a much slower rate than 
the linear growth obtained for usual diffusions. In Harris' original work, a system 
of independent Brownian motions with reflection was considered. These kind of 
systems are known in the physics literature as single-file diffusions (see [6] for a 
recent discussion and further references). 

In this subdiffusive setting, it has been proved that the rescaled position of the 
tagged particle converges, in the sense of finite-dimensional distributions, to a frac- 
tional Brownian motion ([4], [10], [18] for the simple exclusion process; [27] for 
interacting Brownian motions). Recently, a functional central Umit theorem has 
been obtained [26]. This subdiffusive behavior is characteristic of single-file diffu- 
sions; when particles can pass one over the others or in dimension d> I, the tagged 
particle converges to a Brownian motion [20]. For biased particles, the limit is also 
diffusive [21]. 

From the work of Rost and Vares, we know that the displacement xo{t) of the 
tagged particle can be identified as the mass-current through the origin in the in- 
crement process associated to the single-file diffusion. Let {x,(f);/ G Z} be the 
position of the particles at time f > in a single-file diffusion. We assume that 
Xi{t) < x,+i(f). Define rit{i) = x;+i [t) — x,(f) — e, where £ = in the case of par- 
ticles evolving on the real line, and £ = 1 in the case of particles evolving on the 
lattice. The process rjt = {rit{i);i G Z} turns out to have a Markovian evolution 
with a local dynamics. Let Jo{t) be the accumulated current through the bond 
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(—1,0) (of particles in the lattice, and of mass in the continuum). We have the 
identity Jo{t) = xo{t) — xo(0), and therefore the asymptotic behavior of the tagged 
particle is given by the asymptotic behavior of the current /o(0 for the process rjt. 

A second identification, also known from the work of Rost and Vares, allows 
us to obtain the current Jo{t) as a function of the empirical density associated to 
the process T],. At least on a heuristic level, Jo{t) = Li>o{^f(0 ~ ^o(0}- When 
the number of particles is finite, this relation is simply an integrated form of the 
conservation of mass. When the number of particles is infinite, the sum is not 
absolutely summable. We will see that when the process % is in equilibrium, the 
truncated sums Lo<!<n{^<(0 ~ ^o(0} form a Cauchy sequence, and the limiting 
variable is precisely the current Jo{t). 

In this way, the asymptotic behavior of the tagged particle in single-file diffu- 
sions can be obtained in terms of the asymptotic behavior of the empirical process 
associated to the increments of the original process. This approach has been used 
by various authors [23], [27], [6], [14] in order to obtain a central limit theorem for 
a tagged particle in single-file diffusions. 

The idea of relating the position of the tagged particle to the current of parti- 
cles through the origin can also be accomphshed for the original process, without 
considering the increment process. This approach was exploited in great gener- 
ality in [11], where the authors obtain a functional central limit theorem for the 
tagged particle in a system of particles with collisions. Besides the collision rule, 
the evolution is independent. Considering different families of dynamics for the 
motion of one particle, they obtain in the limit any exponent < y < 2 for the 
mean square displacement E[xo{tY] of the tagged particle. The main drawback 
of the approach of [11] is that it does not generalize to systems of particles with 
stronger interaction. 

For fairly general diffusive systems, the so-called hydrodynamic limit of the pro- 
cess r]t is given by a non-linear heat equation of the form dfU = dx{D{u)dxu), where 
D{u) is the bulk diffusion coefficient, which is given by the Green-Kubo formula. 
We say that the process r]t has an hydrodynamic Umit if the rescaled empirical 
process Y^xei, ^tn^{^)^x/n{dx) converges in distribution to a deterministic Umit 
of the form u{t.x)dx, where u{t,x) is the solution of the hydrodynamic equation 
dtU = dx{D{u)dxu). From now on, we focus on lattice systems, so the empirical 
process represents the density of particles in the system. Notice the diffusive time- 
scaling. In particular, if we start the process with a fixed density of particles, at 
a macroscopic level the density of particles does not change. When the invariant 
measures of the process r\t have short-range correlations, the spatial fluctuations 
of the density of particles are of Gaussian nature and they are given, in the macro- 
scopic limit, by ;^(p) #^, where ^ is a white noise in R, p is the density of particles 
and the quantity ;|;;(p) is the static compressibility of the system. By the fluctuation- 
dissipation relation, the density fluctuations around a fixed density p evolves in a 
non-trivial way under a diffusive scaling, and they satisfy in the macroscopic limit 
the infinite-dimensional Omstein-Uhlenbeck equation 



d^ = D{p)mdt + s/D{p)x{p)VdWt, 



(1) 
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where Wt is a space-time white noise. Going back to the representation of the cur- 
rent in terms of the empirical density of particles, we see that putting the Gaussian 
space scaling and the diffusive time scaling together, n^^^^jQ{tn^) should con- 
verge to 'S^tiHo) — '3^o{Ho), where Ho{x) = l(x > 0) is the Heaviside function. 
Since the process is Gaussian, a simple scaling argument allow us to conclude 
that Jo{t) approaches a normal distribution of variance &{p)t^^'^ as t ^ °o, where 

&ip) = x{p)/VD{p). 

Recently ([1], [6]), the following question has been posed. What happens with 
the tagged particle if each particle has its own, different diffusion constant? It turns 
out that [14] under mild conditions, the behavior is the same as before: £'[a:o(?)^] = 
@{py^^, but now the diffusion coefficient D{p) is given by an homogenization 
formula. 

In this article we are interested in systems on which the asymptotic variance of 
the current (and therefore of the tagged particle) is of order t'^ for 7 / 1 /2. From 
a mathematical point of view, the universality of the 7= 1/2 value is pleasant and 
satisfactory, but this is not a good fact from the point of view of modehng. There is 
also experimental and numeric evidence supporting values 7 < 1/2 for the variance 
growth of the tagged particle in some extreme situations. From a physical point of 
view, only values of 7 in (0, 1) are expected. In fact, if at a small time window 
we observe a positive increment on the current Jo{t), this means that the density 
of particles at the right of the origin is larger than the density at left, and therefore 
we expect the current to have negative increments in the near future. This means 
that /()(?) has negatively correlated increments, restricting ourselves to concave 
functions for the variance of /()(?)• 

Notice that the question about current fluctuations can be posed for any one- 
dimensional system, related or not to a single-file diffusion. Therefore, in this 
article we pose the question about the asymptotic behavior of current and density 
fluctuations for general symmetric, one-dimensional particle systems. In particu- 
lar, we propose two classes of models which will allow us to find a central limit 
theorem for 7o(?) in the full range of scales 7 G (0, 1). We recall now that our 
heuristic derivation of the 7 = 1/2 law is very robust, we have only assumed that 
the static fluctuations are normal, that the hydrodynamic limit is diffusive, and that 
the fluctuation-dissipation relation holds. Following the same scheme, we will ob- 
tain a different value of 7 if the hydrodynamic limit of the process r\t holds in a 
non-dijfusive scaling. If this is the case, we say that we are in prensence of an 
anomalous diffusion. In recent works, two symmetric models on which anomalous 
diffusion occurs have been introduced. In [12], a simple exclusion process with de- 
generated bond disorder has been introduced. The simple exclusion process is just 
a system of simple, symmetric random walks in Z, conditioned to never overlap. 
This is probably the most studied example of a single-file diffusion. Bond disorder 
is introduced in the following way. Let {^x'-,x G Z} be a sequence of i.i.d., positive 
random variables. Assume that the common distribution is on the domain of attrac- 
tion of an a-stable law, a G (0, 1). For simplicity, take t,x>l for any x^Z,. This 
does not alter the tail behavior of E,^, which turns out to be the relevant part. We put 
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a wall of size between sites x andx+ 1, meaning that each time a particle tries to 
jump from x to a: + 1 or viceversa (and if the jump is allowed by the exclusion rule), 
the jump is accomplished with probability Due to the heavy tails of the 
dynamics is dramatically slowed down. In fact, the correct time scaling is 
which is always bigger than n^, and the hydrodynamic equation is dtu = dxdwu, 
where W is an a-stable subordinator corresponding to the scaling limit of the walls 
and dw (also denoted by d/dW in the sequel) denotes the inverse of the Stieltjes 
integral with respect toW. As we can see, the randomness of the environment is 
so strong that it survives in the limit, and even the hydrodynamic equation depends 
on the corresponding realization of the environment. This scaling is robust in the 
sense that stronger interaction between particles does not lead to a different time 
scaUng [13]. 

Looking back to the heuristic formula for 7, we obtain 7 = a/(l + a) < 1/2 
for any a G (0, 1), so the current in this model should satisfy E[Jo{t)^] ~ 
for large t. Instead of considering the simple exclusion process with bond disorder, 
in this article we study the zero-range process with degenerated bond disorder. In 
this way we emphasize that the subdiffusive behavior holds regardless of the details 
of the local interaction. For the simple exclusion process with variable diffusion 
coefficient defined in [14], the increment process is precisely a zero-range process 
with bond disorder. Therefore, a central limit theorem for the current through the 
origin leads to a central limit theorem for the tagged particle in this last model, 
which falls into the category of single-file diffusions. 

A second example of a symmetric system with non-diffusive hydrodynamic 
limit is the zero-range process with long jumps [17]. In this system particles in- 
teract between them only when they share positions, and the jump probability of 
the underlying random walk satisfies a power law: p{z) = c/\z\^~^", cc G (0,2). In 
this case the correct scaUng Umit is n", which is always smaller than n^, and the 
hydrodynamic equation is of the form dfU = Aaq){u), where = — (— A)"/^ is the 
fractional Laplacian and <p(-)is a function encoding the interaction between parti- 
cles. In this case, the heuristic formula gives 7=1/0:, which is always bigger than 
1/2. Notice that for a < 1, this formula gives 7 > 1. It turns out that for a < 1, the 
current through the origin is not well defined, since in any time window + 
there is an infinite number of particles crossing from one side of the origin to the 
other in that time window. 

Since the particle jumps are not restricted to nearest-neighbors, it is not possible 
to find a single-file diffusion for which the increment process falls into this class. 
Therefore, for any 7 G (0, 1), we have a model for which the current of particles 
should satisfy E[Jo{t)^] ~ t"^, but for the tagged particle problem we have a model 
for which £'[xo(?)^] ~ only when 7 < 1/2. 

Although the heuristic plan looks simple and it has been accomplished in the 
diffusive case for many examples, anomalous diffusive behavior poses difficulties 
that are absent for diffusive systems. The main technical difficulty relies on estab- 
lishing the fluctuation-dissipation theorem for the process rjt. In the superdiffusive 
case, the main obstacle is that the fractional Laplacian A^ does not leave invariant 
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the Schwartz space ^(M) of test functions, and therefore the classical construc- 
tion of generalized Omstein-Uhlenbeck processes due to HoUey and Stroock [16] 
does not apply. Solutions of (1) can be constructed using the formalism of Gauss- 
ian process for fairly general driving, non-positive operators ^ [3]. However, this 
construction is not suitable for proving convergence when (as it is the case here) it 
is not easy to show that the limiting process has Gaussian distributions regardless 
of initial conditions. Section 2 is basically a recall of [8] and [9], where powerful 
methods have been developed to prove such convergence theorem, very much in 
the spirit of [16]. We include this section with no new results for the reader's con- 
venience, since we are not aware of previous results applying these ideas to lattice 
systems. In particular a notion for uniqueness of equation (1) is stated, and the 
notion of intermediate spaces is introduced. 

In the subdiffusive case, even the definition of the corresponding Omstein- 
Uhlenbeck process poses a challenge. In fact, besides constant functions, there 
are no smooth functions on the domain of the operator = dxdw- Moreover, 
for any two realizations of the subordinator W, the domains of the correspond- 
ing operators have in common only constant functions. The new material 
on this article starts at Section 3. In Section 3 we construct a nuclear Frechet 
space ^ = which will serve as a test space in order to define the general- 
ized Omstein-Uhlenbeck process associated to the operator Our key input 
is a compactness result for weighted-Sobolev spaces, very much in the spirit of 
the definition of the usual weighted-Sobolev spaces in M. Our construction works 
for any increasing, unbounded function W, and could be of independent interest. 
In Section 4 we give detailed definitions of the zero-range process wiht random 
environment and with long jumps. We also state our main results concerning the 
asymptotic behavior of the density and the current of particles. In Section 5 we 
obtain the fluctuations of the density in the superdiffusive case. In Section 6 we 
obtain the fluctuations of the density in the subdiffusive case. A key intermediate 
result is the so-called energy estimate, which roughly says that the space-time fluc- 
tuations of a given function can be estimated by the Dirichlet form associated to the 
underlying random walk. This result holds true for any reversible system, regard- 
less of the super or subdiffusive behavior of the system. The universality of this 
estimate is more evident in Section 8, where we prove fluctuations for the current 
of particles through the origin in both super and subdiffusive cases. We finish this 
article in Section 8 by obtaining the fluctuations of a tagged particle in the simple 
exclusion process with variable diffusion coefficient, as a direct consequence of 
the results in Section 7. We point out that all our results in the subdiffusive case, 
apphes for any process W{x), stochastic or not, such that lim„^±oo = ±00 and 
such that the environment has a version converging almost surely to W{x). 



2. Generalized Ornstein-Uhlenbeck processes 

In this section we give a precise definition to what we mean by a generalized, 
or infinite-dimensional, Ornstein-Uhlenbeck process and we state some conditions 
for uniqueness of such processes. All the material in this section has been taken 
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from [8] and [9]. The interested reader can find a more detailed exposition and 
further apphcations in those articles. 

Remark 1. Throughout this article, we use the denomination "Proposition" for re- 
sult that have been proved elsewhere. We reserve the denomination "Theorem" for 

original results. 

2.1. Preliminary definitions. Let .if^(]R'^) be the Hilbert space of square inte- 
grable functions <p : ^ M. Let ^ : Z)(Jf ) C ^^{W') ^^{R'') be the gen- 
erator of a strongly continuous contraction semigroup {St'J > 0} in ^^{W^). Let 
{| I ■ 1 1„}„ be an increasing family of (not necessarily finite) norms in ^^(M^) with a 
common kernel such that \\(p\\l = f (p{xfdx. Let ^ C ^'^{W^) be the Frechet 
space generated by {|| • ||„}„, that is, the completion of under the metric 

./(<p,VA) = £2-"(||9-VA||„Al). 

n>0 

Let us denote by the topological dual of We can construct in such 
a way that ^ (1 ^^{W^) (1 ^' and such that the inner product (•, •) in ^^(j^rf) 
restricted to #o x #o can be continuously extended to a continuous bilinear form 
(•,•): =^ X M. We assume that the space ^ is nuclear, that is, for any n > 

there exists m> n such that any |j • ||,„-bounded set is a || • ||„-compact set. In 
what follows we will always consider a family of norms || • ||„ for which the set 
{(p G .if^(]R^);||(p||„ < +00} is a Hilbert space under || • ||n, although this point is 
not essential. Our objective is to describe some conditions under which existence 
and uniqueness of solutions can be established for the stochastic equation 

d% = ££*%dt^d%, (2) 

where ^ is a given semimartingale in and ^ is a given, maybe unbounded, 
operator in 5£'^{W^). The canonical example of a nuclear, Frechet space is the 
Schwartz space S^{R^) of test functions in M^. In this case, 

where k denotes a multi-index {k\^... ,kd) with \k\ =ki-\ \-kd and d'^ = Ylidx-- 

The dual of ^(M'^) is the space S^'(W') of tempered distributions and a common 
kernel for each norm 1 1 • 1 1„ is the set "^^{W^) of infinitely differentiable functions 
in M.^ of compact support. 

For a given topological space S", we denote by ^([0,r],<f) the space of cadlag 
trajectories in For simphcity, we consider a finite time interval [0,T]; results 
for [0,00) will follow from standard extension arguments. For a <b,we denote by 
'^°°{a,b) the space of functions in [a,b\ with support contained in {a,b). The 
space ^°°{a,b) is a nuclear Frechet space with respect to the topology of uniform 
convergence on compacts of the function and its derivatives of any order. Notice 
that functions in ''i^°°(a,b) vanish att = a,b. 

Let us denote by #o,r the tensor product ^ ® ^°°{0, T). For a clear exposition 
about tensor products, see [28]. This space is also a Frechet space. We denote by 
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the topological dual of ^o,r- Now we state some technical lemmas which 
will be useful in the identification of solutions of (2). 

Lemma 1. For any trajectory x. in ^([0, T\ , ^'), let x G .^q j. be defined by 



for any Y ^ ^o.r- Then the mapping x. ^ x is a continuous mapping fi-om the 
space &{[0,T],,^') into ^^ j. 

Lemma 2, Let x. be a process in S![ [0, T] , a.s. continuous att = T. Then the 

distributions ofx., x determine each other 

Definition 1. We say that a family {Xf, f G ^} of integrable random variables is 
a linear random functional if Xf is linear and continuous as a function of /. 

Lemma 3. Let Xf be a linear random fiinctional. Then, there exists a unique 
random variable X in ^' such that {X,f) = Xf a.s. for every / G 

2.2. Weak formulation. When the operators JSf and St have the good taste of 
leaving the space ^ invariant, there is a very intuitive notion of solutions for (2). 
This is the case, for example, when ££ = S.. We start doing some formal manip- 
ulations. Here and below, the initial condition ^ will be a random variable with 
values in ^' . Apply (2) to a test function (p, multiply by another test function 
/ G 5,r), integrate the equation over time and perform an integration by 

parts to obtain 



Since we want to capture the initial distribution, we will extend the space ^o,r 
a little bit. Notice that a function / G =^o,r can be thought as a trajectory in ^ 
which vanishes at f = 0, T. Take any 5 > and define ^-sj as the tensor product 
5,r). We define as the set of trajectories in restricted to the 

interval [0, T]. These trajectories always vanish utt = T , and of course ,^o j C 
By linearity and an approximation procedure, we can extend the previous identity 
to test functions \\f G ^t- For ^^f G ^t, this identity reads 



Definition 2. We say that a process ^ in ^([0,7],=^') is a solution of (2) if ^ : 
=^ — > =^ is continuous and (3) holds for any G ^j- 

Let us recall now the variation of parameters method to solve linear evolution 
equations. The equation d'3^t = ■^*'3^t has as a solution the process ^ = S'^&'o, 
which means that ^ is defined via the relation = {'3i),St(p). Notice that 

^ is well defined in ^([0, T],^') if St : ^ ^ ^ is continuous. The variation 
of parameters method suggests to search for a solution to (3) of the form ^ = 
5,*^, where =^ is a semimartingale satisfying = In that case, formally 
satisfies (3) if dJ^t = S%d3ft. Since is not well defined, the usual trick is 





{Wt,(pf(t)+f(t)^(p)dt={%,(pfm+ / {%,(pf'{t))dt. 





(3) 
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to multiply this expression by S*, to obtain a well defined relation for < f < t': 
S*,d% = S*,_^d%. Integrating this expression, we obtain that S*% = S*% + 

So S*_f.d%. 

In terms of test functions, this expression gives 



Jo 



For a test function of the form (p{x)f{t) with (pe^andfe ^°°(0, T), the same 
manipulations yield 

r{m,(pm)dt = {%, r f{t)s,(pdt)+ r dt f {d^s,s,.s(p)m. 

Jo Jo Jo Jo 

The last term can be written as 

I {d%J f{t)S,.s(pdt) = - I {%-(pf{s))-[ f{t)St_,^(pdt)ds 

Jo Js Jo Js 

= - f f'{t)S,.s(pdt)ds. 

Jo Js 

As before, using linearity and an approximation procedure we can extend this 
identity for arbitrary test functions \\f G ^t- 

f {m,¥r)dt = {%, f StWtdt)-[ {2fs,! St-s^\iftdt)ds. (4) 
JO Jo JO Js at 

Definition 3. A process in i^([0, T], J^') is said to be an evolution solution of 
(2) if S( : =^ ^ =^ is continuous, for <p G =^ the curve t ^ (^^, 9) is continuous in 
and (4) is satisfied for any G ^t- 

Notice that by Lemma 2, relation (4) determines the distribution of Now we 
explain on which sense we have a unique solution of (2): 

Proposition 1. Let '3^Qbe a distribution in There exists a unique solution "W. of 
(2) with initial condition if the conditions on ^ and St stated in Definitions 2, 
3 are fulfilled. This solution is given by the evolution solution defined in Definition 
3. 

Uniqueness follows from (4) and Lemma 2. Existence follows from standard 
methods for evolution equations. We call the process the generalized Ornstein- 
Uhlenbeck process of characteristics (^,^) and initial condition The most 
well known example for which this Proposition apphes is the operator ^ = A. 
In that case, it is well known that ^ = S^{W^) is left invariant by and by Sf. 
However, an important example, relevant for our purposes, which falls out of this 
setting is the operator ^ = — (— A)"/^, a G (0,2), i.e., when _Sf is the fractional 
Laplacian in 3^^. The fractional Laplacian is an integral operator which can be 
written in the form 

for a properly chosen constant > 0. Due to the long-range integration, for any 
positive function <p G ^{M.'^), the function ^<p ^ ^{R^): decay at infinity fails. 
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as it is easily shown taking Fourier transforms. Therefore, for the case on which 
jSf = — (— A)"/^ a new interpretation of (2) is needed. This is the content of the 
following section. 

2.3. The intermediate spaces. Throughout this section, we take ^ = ^(M^), 
and therefore ^' = ^'{W'). For p>0 and <p G '^{W'), the space of continuous 
functions in R^, define 

= sup |<pW|(1 +x^)p 

and = {<p G "^{R^y, \\(p\\p,oo < +°°}. Define also 

<^pO = {<p G <^(M^); lim <p(x)(l +x^y = 0}. 

Clearly, ^p, '^pi^ are Banach spaces with respect to the norm || • ||poo. We also 
have the continuous embeddings ^ ^ '^p^^o for p > and "^p^o ^ ^^{R^) for 
p> d/2. Let be the topological dual of ^p,oo. We have the chain of inclusions 

Denote by || • ||-p,oo the dual norm in '^^oo- Let be the set of positive. 
Radon measures fX in such that /(I + x^)^p jx{dx) < +oo. We give to 
the p-vague topology, which is the weakest topology that makes the mappings 
(p<-^ {H,(p) =: J (pdfl continuous for every <p G ^(M'') of compact support and for 
(p{x) = {l+x^)~P as well. Notice that for ^ G^+, ||At||-;,.oo = J{l+x^yPiJ.{dx). 
Notice also that Lebesgue measure is in for p > d/2. Define = — (— A)"/^ 
and let St be the semigroup generated by A^. The following propositions can be 
proved taking the representation of A„ in Fourier space and the representation of 
St in Fourier-Laplace space. 

Proposition 2. For t >0 and d /2 < p < (d + a) /2, St is a bounded linear operator 
from '^p into itself, and also from into itself. 

Proposition 3. Let (p G be such that the limit lim|;f|^t„ <p(x)(l +x^) exists. Then 
t St(p is a continuous trajectory in For any }X G t ^ StjX is a p-vaguely 
continuous trajectory in Moreover, there is a positive constant Cj such that 

for any (p ^'tfp and any <t <T. 

Proposition 4. The mapping <p G =^ i— > A«<p G 'rfpfl is continuous. 

The idea now is to define a generalized version of what we mean by a solution 
of (2). Notice that a priori expresions of the type A^cp) are not well defined, 
since Aa(p does not belong to In the following two definitions, ^ denotes an 

arbitrary Frechet space. 

Definition 4. We say that a process '3^. in ^([0, T] , is a generalized solution of 
the (2) if there exists a Banach space y of functions in such that 

i) ^^r^ ^^{w^) 
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ii) The linear mapping : \% continuous 

iii) jJ{'3^t,^Vt)dt is well defined for any i/A. G {D{^) n y) (g) ^°°(0, T) 

iv) For any Xj/. as in iii), the identity (3) holds a.s. 

Definition 5. A process ^ in ^{[0,T],^') is a generaUzed evolution solution of 
(2) if there exists a Banach space y such that 

i) ^^r^^^iM.'^) 

ii) St : ^ ^ f is continuous for any t G [0, T] 

in) t St(p is a continuous trajectory in ^ for any (p e ^ 

iv) The right-hand side of (4) is well defined for any v^. e ^ (g) ^°°(0, T) 

v) For any l//: as in iv), the identity (4) holds a.s. 

Proposition 5. Under the conditions of Definitions 4, 5, for any initial distribution 
W() in there exists a unique generalized solution of (2), which is given by the 
generalized evolution solution defined in Definition 5. 

Notice that in particular, the stochastic equation d'3^t = ^*a'3^tdt + d% has a 
unique solution for any initial distribution S% in 



3.1. The derivative d/dW . In Section 2 we have discussed existence and unique- 
ness of generalized Ornstein-Uhlenbeck processes defined in dual Frechet spaces. 
Up to this point, the examples discussed so far (corresponding to the fractional 
Laplacians A^) do not require to take an abstract Frechet space, since the only 
space considered was the space ^'(M^) of tempered distributions. In this section 
we give an example where ad-hoc nuclear spaces need to be constructed. 

Let W : M ^ M be a strictly increasing function. By convention, we assume that 
W is cadlag and that W(0) = 0. For two given functions /,g : M — > M, we say that 
g = df/dW if the identity 



holds for any x € M, where the integral is understood as a Stieltjes integral, and 
/o means integration over the interval (0,x]. In the same spirit, we say that g = 



for any x e M. We will denote the operator d/dxd/dW by S^^f. Notice that when 
the function W is differentiable and W\x) ^ for any x (that, W is a diffeomor- 
phism), = dx{W' {x)^^ dx)- We say that a function / is W -differentiable if there 
exists a function g : M ^ M such that g = df /dW . When W is a diffeomorphism, 
this notion reduces to usual differentiabiUty. 



3. Generalized differential operators 




d/dxdf/dW if 
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Now we want to define as an unbounded operator in ^^(M). We say tliat 
a function / G JSf^(R) belongs to D{5^^/f), the domain of ££^ff, if there exists con- 
stants & G M and a function g G ^^(M) such that 



for almost every x G M. In that in this case, / has a cadlag version for which (5) 
holds for any x G M. Notice also that the double integral is well defined for any 
X, since g G ^;j,^(R) and therefore the integrand of the second iterated integral 
is continuous. From now on, we assume that any / G D[^5£yf) is cadlag. If / G 
we define ££^f = g. Evidently, the operator ■ D{^w) ^ ^^(M) is a 
linear, unbounded operator. What it is not so evident, is that is densely defined, 
symmetric and non-positive (remember that Jifw is densely defined if Z)(^v^') is 
dense in =^^(M^)). We state these properties as lemmas. 

Lemma 4. The operator is densely defined. 

Proof. Let C : M ^ [0, 1] be a function such that t;{x)=Qiix(^ (0, 1) and 
C(x) > if X G (0, 1). Fix a<b. Take < 5 < (^7 -a)/2 and define : M ^ R by 

i) for X <a, h{x) = 

ii) for a <x< {a + b)/2. 



whcTcc{b,5) = jlsl^{{b-y)/5)W{dy) 
iv) h{x) = for X > ft. 
By construction, h G D{^w) and when d ^ 0, h converges in ^^(M) to l(x G 
[a,b]), the indicator function of the interval [a,b]. Considering the vectorial space 
generated by the functions h, we conclude that D{^w) is dense in =Sf^(R). □ 

Remark 2. Starting the construction with a function ^ in D{^w) instead of a 
smooth function, we obtain a kernel for (=Sfty)^. A simple iteration of this con- 
struction allows us to prove that in fact the operator {^wT is densely defined for 
any n > 0. 

Lemma 5. The operator S^yv is symmetric and non-positive. 

Proof Take /,g G D{Ji^'w), of compact support. By the previous lemma these 
functions are dense in ^^(R). Then, 




(5) 




whcr^ cia,5)=C^{{y-a)/d)W{dy) 
iii) for {a + b)/2<x<b. 



h{x) =h{{a + b)/2) - c{b, 5)-' / C ((ft - y)/5)W{dy) 
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where we have used the integration by parts formula for Stieltjes integrals and 
the fact that df/dW and dg/dW are well defined and have bounded support. By 

an approximation argument, we conclude the / f^wgdx = / g^-wfdx for any 
G D{J^w), which proves that J^w is symmetric. Taking f = g, v/e see that 
/ fJ/^'wfdx = — J {df /dW)^dW and therefore is non-positive. □ 

3.2. The Sobolev spaces associated to In this section we construct a se- 
quence of nested Hilbert spaces J^n+i ^ --^n • • • ^ Jif^{R) and we define ^ as 
the intersection of such spaces. We will prove that ^ is a nuclear Frechet space 
on which is continuous and therefore we will show that the theory of gener- 
alized Ornstein-Uhlenbeck process of Section 2 applies to ^ and The idea 
is simple. We start recalling the spirit behind the definition of the usual Sobolev 
spaces in R. In finite volume, the Laplacian operator with boundary conditions 
of Dirichlet type satisfies Poincare inequality. From Poincare inequality, we can 
prove that the Laplacian A has a compact resolvent. In particular, the sequence 
of Sobolev spaces associated to the norms = {(p, (— A)"^)^/^ is nuclear, and 
we also know that its intersection contains all the infinitely differentiable functions 
of compact support. However, in infinite volume the operator A does not have a 
compact resolvent, and another definition is needed. We already know that the 
Sobolev norms in infinite volume include a polynomial weight. One way to under- 
stand this polynomial weight is to look at the Helmholtz operator — A + This 
operator has a compact resolvent, and its eigenvectors are the Hermite functions. 
The vectorial space finitely generated by the Hermite functions is the set of func- 
tions of the form p{x)e~'^ , where p{x) is a polynomial. A possible definition for 
tiie Schwartz space ^{R) is tiie following. Take ||<p||„ = (<p, (-A + x^)"<p)i/2, 

2 

^0 = {p{x)e~^ ; p{-)is a polynomial } and define y{M.) as the closure of un- 
der the topology generated by the norms 1 1 • | |n- 

We will proceed in an analogous way for Define ^ = —^w + \W{x)\ 
(here we understand as a multiplication operator). In the following series 

of lemmas, we will construct a nuclear space ^ on which is continuous. We 
assume that W is continuous a.tx = 0. This assumption is not restrictive, since =2V 
remains unchanged by adding a constant to W and W has at most a numerable set 
of discontinuities. Therefore, a simple shift of the origin makes W continuous at 
x = 0. 

We start with a very simple observation. 

Lemma 6. The operators (^nd \W{x)\ have a common kernel ^. 

Proof. It is enough to take ^ as the set of functions constructed in Lemma 4. In 
fact, by construction this set is a kernel for and it is well known that the set of 
functions of compact support is a kernel for any multiplication operator on which 
the multiplication function is locally bounded. □ 

Notice that from this proposition we conclude that the operator ^ is densely 
defined, a fact which is not clear a priori. The same arguments also apply for the 



ANOMALOUS FLUCTUATIONS 



13 



operator s^^, n>\. For <p G D(^") define 

1/2 



{f (<P,^V)}' 



Define now the Sobolev space = as the completion of D{£/") under 

the norm || • ||n. Notice that J^+i C C ■ ■ ■ C ^^(R) and each space is a 
Hilbert space with inner product obtained from 1 1 ■ | |n by polarization. 

Before stating the next lemma, we need a simple definition. We say that a func- 
tion / : M — > M is W-Holder continuous of exponent /3 > if there is a constant 
c> such that \f{x)-f{y)\ < c\W (x) -W {y)\l^ for any x,y e R. 

Lemma 7. Any function (p G D{s^) is W -Holder continuous of exponent 1/2. 

Proof. Notice that 

\m\ = j (^)'^^ + / W{x)\^{xfdx. 



It is enough to see that 

W{y)-<P{x)\^={[^dw)' <\w{y)-W{^ 

<\W\\\\w{y)-w{^)l 



\dW) 



which proves the lemma. □ 

Lemma 8. Let be a sequence in ^^(M) such that ||(p„||i < 1 for every 

n. Then, there are a subsequence n' and a function (p £ Jf^(M) such that ||<p„/ — 
<p| |o — >^ a* n — > oo. In particular, the embedding ^\ ^ .if^(M) is compact. 



Proof By Lemma 7, the sequence is VK-Holder continuous of exponent 1 /2. 

Consider in R the distance d\ff{x.,y') = \W{y) — W{x)\. The metric space {R,dw) is 
not complete, but its completion can be obtained by "splitting in two" each point 
of discontinuity of W. More precisely, let {x,;/ € N} be the set of discontinuities 
of W. Defining Rw = MU {;»:,—},■ and taking dwixi-,y) = \W (xi-) - W {y)\, we 
see that {Rw,dw) is a complete metric space. Now we notice that the functions 
% are 1/2-Holder continuous (in the usual sense) in Rw In particular, we can 
continuously extend (p„ to the whole space Rw. With respect to the usual topology 
of R, this just means that % is cadlag, and that the set of discontinuities of % is 
contained in {x,},. In Rw, the sequence is uniformly equicontinuous, and 

by Arzela-Ascoli theorem, there are a subsequence n' and a function (p such that 
(p„i — > <p, uniformly in compacts. The function (p inherits the Holder-continuity 
from (p„i. Considering cadlag versions of the functions <p„, <p, we see that <p„/ <p, 
uniformly in compacts, also with respect to the usual topology of R. In particular, 
W (p pointwise. By Fatou's lemma, J \W {x)\(p{x)^dx < 1. Since (p is cadlag, 
/ii (p^dx < +°o and we conclude that <p G .if^(M). Since <p„/ (p uniformly in 
compacts, in order to have convergence in ^^(M) as well we only need to have a 
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uniform estimate on the tails of the integrals / (<p„' — (p)^. For any M > 0, 
limsup / {%' — <p) dx< limsup / {(pn' — <p) dx 

fl' ^OO J fl' -rOO J — M 

+ 2 hm sup / {(pl + (p^)dx 

[-M,MY 

^ min{-W(-M),W(M)} /l^(^)l('^"' + ^')-^^ 

^ 4 

- min{-W(-M),W(M)}' 

Since M is arbitrary, we conclude that | |<p„' — <p| |o — as n' — >^ oo. □ 

As a simple consequence of this lemma, we see that the operator ^ has a com- 
pact resolvent. The operator is also symmetric and non-negative. In particular, 
there is an orthonormal basis of ^^(R) formed by eigenvectors of £/. De- 
fine #0 = span{<p,},-, the vectorial space finitely generated by {(pi}. It is clear that 
^0 ^ D{j^") for any n > 0. Since both operators and are local, we can 

take the projection of (pi into a finite box to conclude that (pi is well defined, 
but at this point we can not ensure that ^w^Pi is an element of =Sf^(M). 

Lemma 9. For any function (p G J^i, ^w^> € ^^(M). 

Proof. Let us compute (^^<p, <p) for <p e D{s^): 

= = w^wfpW^ + \\\w\^''^(p\? + '^{--^w(p,\W\(p). 

Let us compute the crossed term (— =Sfty<jp, \W\(p) with care. By Remark 2, there 
is a kernel for s^^ composed only of functions with compact support. Take a 
function <p G D{j^^) of compact support. In that case, 

J ' '\dW I J IdW ' ' 

where we have performed an integration by parts in the second Une and we have 
used Leibniz's rule d{fg) = fdg + gdf in the third Une. Therefore, 



(^>,<p)=ii^w<pii§+iiiw|V>i|2+2|||wr/2j<p/jwii^-(<p(o-)2+<p(o)2). 
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In particular, we have proved that 

I \^w(p\\l < (=2/'<P, <P> + <p(0-)2 + <p(0)2 

< \\(p\\l + (p{0-f + (p{0)^. 

Therefore, our task is now to bound 9(0) and 9(0—) in terms of 1 19| I2. Since W 
is continuous atx = 0, by Lemma 7 <p(0— ) = <p(0). If <p(0) = 0, there is nothing to 
prove. Assume without loss of generality that (p (0) > 0. Take a strictly increasing, 
continuous function T : [0,oo) ^ [0,oo) such that r(0) = and \W{x)\ < r{\x\)^ for 
anyxeR. Since \(p{x) - (p{o)\ < ||9||i|W(;c)|i/2, taking § = r-\(p{0)/2\\(p\\i) 
we have that <p{x) > <p(0)/2 for x G [-5, 5]. Therefore, 



2, 9(0) 



x) dx 



) ^ 




Define &{x) = x/2^/^^^x/2). Notice that : [0,°°) [0,°°) is continuous 
and strictly increasing. Therefore, (p{0) < ||9||i0~^(||(p||o/||(p||i). Noticing that 
l!<Pl!{)< 1 1 9 1 1 1, we have proved that 9(0) < ||9||o0^^(l). Putting this estimate into 
(6), we conclude that there exists a constant c = c{W) such that 

11^1^9110 <c(W)||9||2 

for any function of compact support. Taking suitable approximations, the lemma 
is proved for any function 9 G J^. □ 

Notice that the previous lemma is saying something more about In fact, we 
can conclude that is a continuous operator from to =Sf^(M). Define in ^0 
the metric 

d{f,g) = l^2-"\\f-g\UAl 

n>0 

and let ^ be the closure of ^0 under this metric. By definition, ^ C Jif„ for any 
n and by Lemma 8 the space ^ is a nuclear Frechet space. 

Lemma 10. The operator ^yi/ maps continuously ^ into itself. 

Proof. Repeating the arguments of Lemma 9, we see that JSfiy maps continuously 
into for any n > 0. □ 

Now we are finally in position to prove the main result of this section. 

Theorem 1. Let 2ft be a semimartingale in , the topological dual of ^ with 
respect to J^^ (R), and let "S^obea distribution in ^' . There exists a unique solution 
'3^t of the equation 

d% = ^^'3/tdt^d% 

with initial condition 

Proof We just need to check if the conditions of Proposition 1 are fulfilled by 
iri Since the operator is bounded in the semigroup St = e'-^^ is 
well defined by taking the corresponding power series, for example. For a given 
9 G the trajectories 5^9 are continuous in ^ by the semigroup property and 
the boundedness of JSf^y, so Proposition 1 applies in this situation. □ 
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4. The zero-range process 

Let p = {p{x,y);x,y ^l/}be the jump rate of a continuous-time random walk 
on U^. Let g : No — [0)°°) be such that g(0) = and g{n) < n for any n G Nq. 
The zero-range process with jump rate p and interaction rate g{-) is a system of 
random walks on with the following dynamics. For each ordered pair of sites 
(x,_y) G Z'^ X Z'^, after an exponential time of rate g{n)p{x,y), where n is the number 
of particles on site x, one particle jumps from site x to site y. Notice that the 
condition g{0) = ensures that there is at least one particle at site x when the 
jump occurs. After a jump, a new exponential time starts afresh. This happens 
independently for each pair of sites {x,y). 

Since p is the jump rate of a continuous-time random walk on Z'^, we have 
LyGZ'' Pi^^y) < +°° for X £ Z"^. Without loss of generality, we assume that 
p{x,x) = for any x G Z^. Since we are interested in symmetric systems, we also 
assume that p{x,y) = p{y,x) for any x,y G Z'^. Notice that the dynamics described 
so far is conservative, that is, particles are not created neither annihilated. When 
the initial number of particles is finite, this dynamics corresponds to a continuous- 
time Markov chain. Let H = Nq be the space of configurations of particles for 
this process. We denote the state of the process at time t by rjt, and we denote 
by rit{x) the number of particles at site x at time t. When the initial number of 
particles is infinite, some conditions are needed in order to prevent explosions on 
the system, that is, the arrival of an infinite number of particles to some site in finite 
time. For our purposes, the Lipschitz condition sup„ \g{n+l) — g{n)\ < +°o on the 
interaction rates and sup^£'[T7 (x)] < +oo on the initial (maybe random) distribution 
of particles will be enough (see [2], [24] for a more detailed discussion). 

The process J]f described above is a Markov process generated by the operator 
L given by 

where / : H ^ M is a local function (that is, it depends on T] (x) for a finite number 
of sites X G Z'^), 7]^'^ is given by 

{T](x)-1, Z=X 
r]{y) + h z = y 
r]{z), zj^x,y 

and / satisfies a Lipschitz condition: 

sup|/(T]"'^)-/(r7)|<+oo. 

4.L Invariant measures. Under the symmetry condition ;7(x,y) =;7(y,x) for any 
x,y, the process % has a family of invariant measures of product form. For <p > 0, 
let us define 
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where g{n) ! = g(l) • • • g{n) for « > 1 and g{0) ! = 1. Define 

(Pc = liminf </g{n)\. 

n— >oo 

The power series defining Z((p) converges for (p < (pc and diverges for (p > (pc. 
We assume that (pc > 0. Tliis is the case if, for example, liminf„g(«) > 0. For 
<p < (pc, let Vfp be the product measure on D. defined by 

v,(r](xO=ni,...,n(x,)=nO=n^^- 

The measure v,p satisfies the detailed balance condition, and therefore it is left 

invariant by the evolution of rjt and it is also reversible under rjt. Let us de- 
fine the density of particles p{(p) = J ri{x)V(p. A simple computation shows that 
p{(p) = (pZ'{(p)/Z{(p). Since (p p{(p) is strictly increasing and smooth, Pc = 
Um,p^,p^p(<p) exists and p : [0,(pc) [0,Pc) is a diffeomorphism. 

Since the number of particles is conserved by the dynamics, it is more natural 
to parametrize the family of invariant measures {V(p;(p G [0, (p^)} by the density of 
particles. Denote by <p(p) the inverse of p(<p). We define Vp = v^(p). 

The family {Vp;p G [0,Pc)} of invariant measures is in addition ergodic if the 
transition rate p is irreducible and g{n) > for n > Q. From now on, we assume 
that in fact the measures {Vp} are ergodic under the evolution of rjt. 

4.2. Density fluctuations. Let us fix p G (0,Pc) and let us take Vp as the initial 
distribution for Letx(?) be the continuous-time random walk associated to p{-). 
Assume that x{t) satisfies an invariance principle in the following sense: there are 
a scaling sequence {a{n);n G N} and a Markov process ^(?) such that 

lim -x(ta(n)) = 3^{t), 
n— >°° n 

in distribution with respect to the Skorohod topology on the space of cadlag paths 

^([0,oo),M'^). Let us denote by ^ the generator of the process S^{t). Assume that 
there are a nuclear Frechet space ^ and a Banach space y for which the conditions 
in Proposition 5 are satisfied. Let '^c(lR^) be the space of continuous functions in 
of compact support. Since the measure Vp is of product form, for any function 
G G 'ic{^'^), 

1™ ^ L nt{x)G{x/n)=p / G{x)dx, 

in probability with respect to the distribution of r\t. It is natural, therefore, to define 
the density fluctuation field as the process in ^([0,°°),^') given by 

^"(^) = 42 I ij\tain){x)-p)G{xln) 

for any G G Since the measure Vp is of product form and ^ (x) has a bounded 
second moment with respect to Vp , for each fixed time f > 0, the field converges 
to a Gaussian field of mean and covariance matrix ;c(p)5(j — x), where x(p) = 
Var(7](x); Vp) and b(x) is the Dirac mass at x. In other words, converges to 
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a spatial white noise. The fluctuation-dissipation principle states (at least in the 
diffusive setting) that there exists a constant D{p), called the diffusivity, such that 
converges in distribution to the solution ^ of the Omstein-Uhlenbeck equation 

d^,=D{p)^*^,dt + d2ft 

with initial condition equal to the spatial white noise described above, where 
^is an ^'-valued martingale of quadratic variation (^) characterized by 

{m){G) = mG))=2D{p)z{p)t{G-^G). 

In the context of the zero-range process, D{p) = (p'{p) and ;^(p) = (p{p)/(p'{p). 
Our aim is to give a rigorous proof of this principle in two situations: for a superdif- 
fusive system on which a{n) grows faster than and for a subdiffusive system on 
which a{n) grows slower than the diffusive scaUng n^. 

4.3. The superdiffusive case. In order to simplify the exposition and to concen- 
trate in the technical problems posed by anomalous diffusive behavior of the asso- 
ciated random walks, from now on we assume that g{n) = I for any n > 0. Our 
first choice for the jump rates is p{x,y) = q{y — x), where ^ : \ {0} {0,o°) is 
a continuous function satisfying q{Xx) = for any x G M''\ {0}, any 

A 7^ and some a G (0,2). The condition a > is assumed in order to have 
Y,yP{x,y) < +°°- In the other hand, the condition a < 2 is assumed in order to 
have a non-Brownian scaling limit for the associated random walk. This choice 
of jump rates corresponds to the zero-range process with long jumps introduced in 
[17]. It is well known that the associated random walk x{t) satisfies an invariance 
principle: 

limx(m«) = jrm, 

n— >oo 

where ^{t) is a Levy process generated by the operator ^ whose action over 
smooth functions F of compact support is given by 

^F{x) = J q{y){F{x + y)+F{x-y)-2F{x))dy. 

Since ^(O is continuous and positive in the unit sphere, we see that =Sf is bounded 
with respect to the fractional Laplacian — (— A)"/^. In particular, ^ also satisfies 
the conditions of Proposition 5, and therefore there exists a unique generalized 
solution of the Omstein-Uhlenbeck equation 

d^t(p'{p)^*^tdt + d3ft (7) 

with initial distribution equal to a mean-zero spatial white noise of covariance 
matrix ;t(p)5(y — x) and driving martingale 3ft satisfying 

(^(G)) =2?9(p) JJ dxdyq{y){G{y)-G{x))\ 

A simple integration by parts shows that this last integral is equal to the energy 
form (f(G,G) =: — j G{x)^G{x)dx. Denote by T]/' the zero-range process r\t„a 
associated to this choice of jump rates and with initial distribution 7]q given by Vp . 
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We adopt the notation P" for the distribution of the process rj" in ^([0, T],D.) and 
E" for the expectation with respect to P". 

Theorem 2. The fluctuation field defined by 

converges to the Ornstein-Uhlenbeck process solving (7), in distribution with re- 
spect to the Skorohod topology in &{[0,T],^'). 

4.4. The subdiffusive case. In this section we define a process evolving in the 
one-dimensional lattice. A positive random variable ^ is said to be a-stable {a G 
(0, 1)), if it satisfies 

log£[e-^?] = -c|A|« 
for some constant c > 0. Let {^x;x e Z} he a sequence of i.i.d. random variables 
with common distribution ^, where ^ is an a-stable law. 

For a given realization of {^x}x, we define p{x,x+ 1) = p(x+ l,x) = \ and 
p{x,y) = if |x — y| 7^ 1. It is well known that the sequence {^x}x satisfies an 
invariance principle: 

1 H 

in distribution with respect to the Skorohod topology of ^(M,M), where W{x) is 
a double-sided, a-stable subordinator with W{0) = 0. Here we denote by [nx] 
the integer part of nx. The following properties of the process W{t) hold a.s. 
The function x ^W{x) is cadlag, increasing and of pure-jump type. We have 
Mm„^±^W{x) = ±oo. The function W{x) is continuous at 0, that is, is not a point 
of jump ofW{x). 

Consider the random walkx(/) associated to the jump rate p{-, •) defined above. 
The process x{t) satisfies an invariance principle [19], [12]: 

Proposition 6. There exists a self-similar process ^ (?) of continuous paths, with 
self-similarity index a/(I -|- a) such that 

lim -x(mi+i/«) = ^(t) 
in distribution with respect to the Skorohod topology o/^([0,oo),M). 

We point out that this is an annealed result: the convergence result is true only 
when averaged over the environment £, . We will not give an exphcit construction of 
the process here; for a detailed description of ^{t), see [19] or [12]. We just 
point out that ^{t) is not Markovian, due to the averaging over the environment. 
Notice the subdiffusive scaling 1 + 1 /a > 2. 

This proposition has been obtained by considering a different environment, dis- 
tributed as £, , but changing with the scaling. A simple computation shows that the 
distribution of IV(1) is equal to A scaling argument shows that, in fact, any in- 
crement of the form £^^/" (W{x + £) — W{x)) is distributed according to ^. More- 
over, if the intervals (xo,xi], (yojji] are disjoints, the variables — W{xo), 
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W{yi) — W{yo) are independent. In view of these observations, for each n we de- 
fine 

^;: = n'/-{W{ix+\)/n)-W{x/n)}. 

Observe that for any fixed n, the sequence ^" = {^"■,x e Z} has the same distri- 
bution of £, , that is, it is a sequence of i.i.d. random variables with common distri- 
bution ^. In the other hand, by construction the normalized sums of the triangular 
array {^";x G Z,n > 1} converge a.s. to the process W{-). As above, we define the 
jump rate Pn{-,-) by taking pn{x,x+l) = pn{x + l,x) = {^")~K and /7„(x,y) =Oif 
\x — y\ / 1. Denote by the random walk associated to p„{-, •). The sequence 
of random walks {xn{t)}„ satisfies the following invariance principle: 

Proposition 7. Fix a realization of W{-). Then, 

lim -;c„(mi+i/«) = ^w(t) 

in distribution with respect to the Skorohod topology of ^{[0,o°),'K), where ^w{t) 
is the Markov process generated by the operator .5^ = d/dxd/dW. 

The main difference between this proposition and the previous one is that now 
the statement of this proposition is a quenched statement, since the random envi- 
ronment is first fixed, and then an invariance principle is claimed for that particular 
choice of the environment. We point out that the process ^w{t) is Markovian, but 
not strongly Markovian. For this and other properties of the process ■^(?), see 
[12]. 

Notice that almost any realization of W(-) satisfies the assumptions made in 
Section 3. Therefore, Theorem 1 holds for the generator and there is a unique 
solution of the Omstein-Uhlenbeck equation 

d'3^^ = (p'{p)^w^^dt + (p{pYl^d^^ (8) 

with initial distribution ^ equal to a mean-zero spatial white noise of covariance 
matrix X{p)S{y — x) and driven by the martingale satisfying 

Let 7]"'^ be the rescaled zero-range process I7f„i+i/a associated to the transition 
rate Pn{-,-)- We call this process the zero-range process with environment W. 

Theorem 3. Fix a realization ofW and a density p G (0,Pc). Consider the zero- 
range process T)"'^ with initial distribution Vp. The fluctuation field ^"'^ defined 
by 

^"'^(G) = ^£(r?r(x)-p) 

converges to the Omstein-Uhlenbeck process '3^^ defined above, in distribution 
with respect to the Skorohod topology in ^([0,°°),=^'). 
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Notice that an annealed result for the process T]"' imphes an analogous result 
for the process t]". In order to state a fluctuation result for the empirical density 
associated to r\f, we need first some definitions. The main difficulty is that the 
space where '3^^ is defined depends on W . Remember that for any W, any / > 
and any G G '^c(IK) the random variable [G] is well defined. In fact, the finite- 
dimensional distributions [G) , . . . [G)) are also well defined. What is 
not well defined, is ^^(G) as a process with cadlag paths. By Kolmogorov's 
consistency theorem, we can consider 'S/'^ as a process defined in ^(M)[°'°°\ 
where ^(M) is the space of Radon measures in M. This space is very irregular, 
but it has the advantage that it does not depend on W . Let P*^ be the distribution 
of in J^{^f'-'"\ We denote by ^™ the process defined in ^(]R)["^°°) and 
having law PCS P*^. The following result is an immediate consequence of Theorem 
3: 

Corollary 1. Fix p G (0,Pc). Le? T]^ he the zero-range process with jump rate 
p^ {■■,■) and initial distribution Vp. The fluctuation field defined by 

^t^'iG) = i= £ {nl,,iM-p)G{x/n) 

for G G "^c^^) converges in the sense of finite-dimensional distributions to 

As the invariance principle for x{t) discussed before, this statement is an an- 
nealed result, which holds only after averaging over the environment. As we dis- 
cussed above, the nuclear space ^ constructed in Section 3 depends on the realiza- 
tion of W. In particular, we do not expect a functional central limit theorem to hold 
for since there is no good set of test functions G that works for any choice of 
the environment. Notice that we have constructed the process only to obtain 
this corollary, which does not involve the test space . 

4.5. Fluctuations of the current of particles. In this section we restrict ourselves 
to dimension d= \ and we consider either the zero-range process with long jumps 
for 1 < a < 2 or the zero-range process with environment W for < a < 1. Given 
two different sites x,y G Z, we define Jx,y{t) as the total current of particles between 
X and y up to time /. That is, Jx,y{t) counts the number of jumps from xioy minus 
the number of jumps from y to x, up to time ? . It is not difficult to see that Jx,y (?) is a 
well defined Poisson process, regardless of the dimension or the parameters of the 
zero-range process. With our definitions, Jx,y{t) = —Jy,x{t), and for the zero-range 
process with environment W, Jx,y{t) = unless x,y are neighbors. The conservation 
of the number of particles can be stated as a continuity equation for T\t{x): 

nt{x)-nQ{x) = -Y,Jx,y{t). 

For X G Z, we define Jxit) as the total current of particles through the bond 
(x,x+ 1), that is, 

y<x 
z>x 
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The inhomogeneous Poisson process Jx,yit) is a birth and death process with 
birth rate p{x,y)g{rit (x)) and death rate p{x,y)g{rit (y)), and therefore the birth and 
death rates of the process Jx{t) are respectively 

£ P{y, z)g{r]t (y)), £ p{z,y)g{r]t {z}). 

y<x y<x 
z>x z>x 

For the zero-range process in random environment W , these two sums are equal 
to pn{x,x+ l)g{r\t{x)) and p{x,x+ \)g{r]t{x+ 1)) respectively. For the zero-range 
process, these sums are convergent for "reasonable" configurations of particles (for 
example, a.s. with respect to Vp) precisely for a > 1. When the number of particles 
is finite, the continuity equation tells us that 

Jx{t) = Y,{n,{y)-no{y)), (9) 

y>x 

that is, the current through the bond 1) up to time t is equal to the number of 

particles to the right of x at time t, minus the number of particles to the right of x 
at time 0. Although equation (9) does not make sense for the process starting from 
the measure Vp, let us assume for a moment that (9) holds in that case. We will 
have 

^Htn^) = i= £ (r7,(,) - r^oiy)) - (^o"(3') - p) 

where Hq{x) = 1(x > 0) and 7 = 1/a for the zero-range process with long jumps 
and 7 = 1 + 1/a in the case of the zero-range process with environment W. Of 
course, ^"(i/o) is not well defined, but at least heuristically we can guess that 
n-^/Vo(m^) converges as « ^ 00 to %{H()) - %{Hq). 

The key idea introduced by Rost and Vares [27] is that the difference '3^" {Ho) — 
^q{Hq) is well defined and it is well behaved with respect to the limit in n. The 
asymptotic Umits for the current are the following: 

Theorem 4. 

i) For the zero-range process with long jumps starting from the initial distri- 
bution Vp, 

lim -7o(m2") = -^(0 
n— »°° n 

in the sense of convergence of finite-dimensional distributions, where ^ it) 
is fractional Brownian motion of Hurst exponent 1 /2a. 

ii) For the zero-range process with environment W and starting from the initial 
distribution Vp, 

hm l7o(m2+2/«) = ^^(t) 

n—nx> fi 

in the sense of convergence of finite-dimensional distributions, where ^{t) 
is a Gaussian process with mean zero and covariance matrix depending on 
W. 



ANOMALOUS FLUCTUATIONS 



23 



iii) Let us denote by Jq (t) the current associated to the zero-range process 
with random environment B,. Under the conditions ofii), 

lim -7o(f«^+^/") = ^{t) 

n— »oo fi 

with respect to the annealed distribution IP" and in the sense of finite- 
dimensional distributions, where ^{t) is a fractional Brownian motion of 
Hurst exponent a/(2 + 2a). 

5. Density fluctuations: the superdiffusive case 

In this section we prove Theorem 2. The proof follows the usual approach to 
convergence of stochastic processes, which consists of two steps. First we prove 
that the sequence of distributions of {^."}n in ^([0,7],^') is tight, from which 
we conclude that the sequence has at least one accumulation point '3(, defined as a 
cadlag process on After that, we prove that the process "Sds a solution of the 
Omstein-Uhlenbeck equation (7). By Proposition 5, there is at most one of such 
solutions. We finally conclude that the sequence is relatively compact and 

has only one limit point, and therefore it is convergent. 

5.1. Tightness of In this section we prove tightness of {'3^"}„. To begin, 

we recall Mitoma's criterion [25]: 

Proposition 8. Let H be a separable Hilbert space, and let ^ C. H be a nu- 
clear Frechet space. Let ^' be the topological dual of ^ with respect to H, and 

let {=S"}„ be a sequence of distributions in &{[0,T],,^'). For a given function 
Ge^, let ^"^^ be the distribution in ^{[0,T],R) defined by ^"■<^(>'(-) € A) = 
(■){G) G A). The sequence {=S"}„ is tight if and only if{^"'^}n is tight for 
any G G J^. 

Applied to the sequence {'3("}n, this theorem says that is tight if and 

only if is tight for any G G y{W^). By Dynkin's formula and after 

some manipulations, we see that 

^,"{G) = %"{G)+ f4j2 I {gK{x))-9{p)}^nG{x/n)ds + ^,"{G), 

(10) 

where J^p{G) is a martingale of quadratic variation 

(^;(G))=2 i: q{}^)iG{y/n)-G{x/n)fg{r]^{x))ds. 

Here we use the notation ^„G for the discrete approximation of ^G given by 
Jf„G(x) = 4 £ q{y/n) (G(x + y/n) + G{x - y/n) - 2G{x)) . 

Notice that J^,G(x) is just a Riemann sum for ^G(x), and due to the symmetry 
of the kernel q{-) and the smoothness of G, the sum converges to ^G^x), uniformly 
inx. 
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In order to prove tightness for the sequence {^?("(G)}„, it is enough to prove 
tightness for {%"{G)}n, {^"{G)}„ and the integral term in (10). The easiest 
one is the initial condition: we have already observed that ^"(G) converges to a 
normal random variable of mean zero and variance x{P) J Gixf^dx. For the other 
two terms, we use Aldous' criterion which now we state. 

Proposition 9 (Aldous' criterion). A sequence of distributions {^P"} in the path 
space ^([0,r],R) is tight if: 

i) For any t G [0,r] the sequence {^"} of distributions in R defined by 
^^(A) = ^"(yit) eA) is tight, 

ii) For any £ > 0, 

limlimsup sup [|j(t + 0) - j(t)| > e] = 0, 

S>Q n— >oo TgTj, 
0<S 

whereTi is the set of stopping times bounded by T and y{z + 6) =y{T) ifz + d > 
T. 

Now we prove tightness of the martingale term. By the optional sampUng theo- 
rem, we have 

P« [| - ^."(G) I > e] < ^E" [< A%(G)> - <^."(G)>] 

20 

<—^„iG,G), 

where in the last inequality we have used the fact that g(-) is bounded, and we have 
denoted by <?„(G,G) the Dirichlet form associated to J5f„: 

4(G,G) = 4 I ^(^)(G(y/n)-G(Vn))'. 

In particular, the martingale J^p{G) satisfies the conditions of Aldous' criterion. 
The integral term can be treated in a similar way: 



^ I {g{n^{x))-^{p))^„G{x/n)ds) 



< 



<0 rE«[(42 I (g(i7;W)-9(p))^nG(x/n)) 
< ^ I Var(g(r](x)),Vp)(^„G(../n))' 



ds 



n 

XI 



<C6T 

for some constant C which depends only on G and p. Therefore, we conclude that 
the sequence {^^"}„ is tight. Notice that as a by-product of these computations, the 
sequence of ^'-valued martingales {^"}n is also tight, as well as the sequence of 
integrals. 
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5.2. The Boltzmann-Gibbs principle. Let "3^. be a limit point of the sequence 
{'3("}n- We want to prove that ^ is a generalized solution of (7). It is enough to 
prove the following result. 

Theorem 5. For any function Ge 
i) The integral 



T 





is well defined as a random variable in ^' . 
ii) The process 

J(,{G) = ^r{G) - %{G) - (p'ip) f^,{^G)ds 

Jo 

is a martingale of quadratic variation {^t{G)) = 2t(p{p)£'{G,G). 

We will assume this theorem and we will finish the proof of Theorem 2. We 
have already seen that for any fixed time t, the =^-valued random variable is 
a Gaussian field of mean zero and covariance matrix ;f(p)5(y — x). Remember 
that {^"}n is tight. Taking an adequate subsequence, we can assume that there is 
an ^'-valued process for which (G) is as above. Since all the projections 
of are martingales, we conclude that is an =^'-valued martingale. Taking 
piecewise linear approximations, it is not difficult to see that for any Xj/t G 

^tiVt) - %{¥o) - J^^ ^s[^^Ws + ^w) ds 

is a martingale of quadratic variation 2x{p) Jo'^{VsjVs)ds. This martingale is 
simply equal to the stochastic integral 

^ d^s{Vs)ds = ^tiVt) - J^s{^^W^ds. 
For t = T v/e have 

%iYo) + ^t[j^Vt+^Vt)dt = ^t[^^¥t)dt, 

which proves that ^ is a generalized solution of (7). 

Before entering into the proof of Theorem 5, we state a proposition which is 
known in the literature as the Boltzmann-Gibbs principle, which was introduced 
by Rost [5]: 

Proposition 10 (Boltzmann-Gibbs principle). Let H -.W^ be a continuous 
function of compact support. Then for any t G [0, T], 



lim E" 



^'4j2L {gK{x))-<pip)-<p'{p)mx)-p))H{x/n)dsy 



= 0. 



The idea behind this proposition is the following. Fluctuations of non conserved 
quantities are faster than fluctuations of conserved quantities (the number of parti- 
cles in our case). Therefore, when we look at the right scaling and when averaged 
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over time, only the projection of the fluctuations over the conserved field (the den- 
sity field in our case) are seen. For the case of the zero-range process with long 
jumps, the proof of this proposition follows the proof for the diffusive case with 
slight modifications, so we omit it and we refer to Chapter 11 . 1 of [22] for a proof. 
We just point out an ingredient of the proof that can be easily overlooked, which 
is the ellipticity condition inf{^(z); \z\ = 1} > 0. This condition allows to compare 
the dynamics with the one associated to the fractional Laplacian A^. We will give 
a more detailed discussion about the Boltzmann-Gibbs principle in Section 6. 

Proof of Theorem 5. We start by proving /). The main point is that J^G is not in 
^ and by tightness we only know that '3^t{H) is well defined for H e Our 
first step is to prove that the operator Jq '3^s{^G)ds is continuous as a hnear 
functional in o (remember the condition d/2 < p < {d + Cc)/2). The operator 
^ is continuous from ^ to 'tfpfi. But 



E"[(^"(7/))^] = 1 £ Hix/nfxiP) < ^Mffl^ (l + ix/nf)-" 



Therefore, for a sequence {^;}/ of functions in ^ which is of Cauchy in o> 
the sequence {/q '3(^"(H)ds} is a Cauchy sequence in ^^(P"), uniformly in n. We 
conclude that the random functional /q &'"{H)ds extends to a random, continuous 
functional in '^pfi. Up to here, we have proved that the family of random variables 
{/o '3^s{^G)ds',G G ,9^} is linear and continuous. By Lemma 3, there is a unique 
^'-valued random functional associated to this family, which proves /). 

Now we continue proving //). Using the Boltzmann-Gibbs principle into (10), 
we see that 



plus a rest that vanishes Remember that our proof of tightness 

actually proved that all the terms above are tight. Taking further subsequences if 
necessary, we obtain that 



for some ,_^'-valued process Notice the convergence of the integral term does 
not follows directly from the convergence of but from /). The convergence 
in distribution of the .^'-valued processes {^"}n does not allow to conclude that 
J^. is an ^'-valued martingale. Notice the simple bound sup„E''[^y (G)^] < +0° 
for any G G We claim that this bound is enough to conclude that is a 
martingale. In fact, from this bound we conclude that the sequence {^j{G)}n 
has a subsequence converging to J^t{G) with respect to the weak topology of 




<C{p,p)\\H\\p^^. 
This computation, when applied to the integral /q '3(l^{H)ds gives 
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^^(P"). Therefore, for a measurable set U with respect to the canonical a- 
algebra 5f, E"[^;'(G)1[/] converges to £[^f(G)lj/]. Since ^"(G) is a mar- 
tingale, E"[.^^(G)lj/] = W[Ji^p{G)lu]. And taking a further subsequence if nec- 
essary, this last term converges to E[^t{G)\u\, which proves that ^.{G) is a 
martingale for any G^^ . 

In order to finish the proof, we only need to obtain the quadratic variation of 
^t{G). A simple application of Tchebyshev's inequality proves that (^"(G)) 
converges in probability to 2t(p{p )S'(G,G) asn^o°. Remember the definition of 
quadratic variation. We need to prove that ^f(G)^ — 2(p{p)tS'{G,G) is a mar- 
tingale. The same argument we used above applies now if we can show that 
sup„E"[^f (G)"^] < -Foo and sup„E"[(^f (G))^] < -^oo. Both bounds follows eas- 
ily from the explicit form of (^," (G)) and (10). □ 



In this Section we prove Theorem 3. The scheme of the proof is the same we 
followed for the superdiffusive case. That is, we first prove that the fluctuation field 
{'&(""^}n is tight and then we prove that any limit point of {'3(""^}n is a solution 
of the generalized Omstein-Uhlenbeck equation (8). With respect to the superdif- 
fusive case, we point out two differences. First, the operator is continuous in 
the ad-hoc nuclear space and therefore the conditions required for uniqueness 
of solutions of (8) are simpler. In the other hand, the environment is degenerated, 
in the sense that the trajectory W has a dense set of jumps and therefore there are 
rates <^ " arbitrarily small at any scale and any macroscopic region. This fact makes 
the system less ergodic, and therefore the proof of the Boltzmann-Gibbs principle 
will be more involved. 

From now on we take a fixed realization of the environment W, so we will omit 
in the notation the dependence in W of different quantities we will define. In order 
to profit from the concepts already introduced in the superdiffusive case, we denote 
by the same symbols analogous objects, like the generator jSf = d/dxd/dW , the 
rescaled process T]", etc. The reader should use the corresponding subdiffusive 
definitions of the various quantities of interest, instead of the ones used in the 
previous section. 

6.1. Tightness of {^"}n. Before entering into the proofs, we need some defini- 
tions. For each n > and each function G : R — ^ M we define 



The operator ^„ is a discrete approximation of ^ = d/dxd/dW and corre- 
sponds to the generator of the process n~^x„(?n^"'"^/"), where x„(?) is the random 
walk associated to pn- It will be useful as well to define the energy form 



6. Density fluctuations: the subdiffusive case 



^nG{x/n) = n 



G{{x+l)/n) -G{x/n) G((x/«) - G((x - \)/n) 1 
W{{x+\)/n) - W{x/n) ~ W {x/n) -W{{x-l) /n) J' 



<f„(G,G)=2^ 



G{{x+\)/n) -G{x/n) 
w{{x+\)ln)-W{x/n 




{w{{x+\)/n) -W{x/n)^ 
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In the definition of S'„(G,G) we have not cancelled out the differences W({x + 
l)/n) — W{x/n) in order to stress that S'„{G,G) is nothing but a Stieltjes sum 
for J{dG/dW)^dW . In particular, for G G =^ the energy S'„{G,G) converges to 

(G,G) = J{dG/dWfdW as « ^ oo. 

The proof of tightness for is very similar to the proof in Section 5.1. By 

Mitoma's criterion, it is enough to prove tightness for {'3("{G)}„ for any G G 
By Dynkin's formula, 

W,"{G) = %"{G)+ 52 (g(77;(x))-<p(p))jSf„G(x/n)J. + ^;(G), (11) 

JO n I rjj 



where ^"(G) is a martingale of quadratic variation 



Tightness of the martingales {^"(G)}„ follows from the deterministic bound 
(^"(G)) < 2tSn{G, G). Tightness of the integral term follows after approximating 
^nG{x/n) by £^'G{x/n) and using that ^G is regular and square integrable. Since 
converges to a white noise of covariance matrix x{P)^iy ~ we conclude 
that {^^"}„ is tight. 

6.2. An estimate of energy type. In this section we will assume the Boltzmann- 
Gibbs principle. We will prove that the limit points of are thus concentrated 

on solutions of (8). A key energy estimate (see Theorem 6) will be the main tool 
allowing to prove the characterization of limit points of {^"}n, explaining the 
name of the section. 

In the statement of the Boltzmann-Gibbs principle (Proposition 10), we ask the 
test function H to be continuous and of compact support. This point is crucial in the 
proof of Proposition 10, since it allows to perform a space average that in the end 
reduces the proof to a version of the ergodic theorem. First of all, it is not clear if 
there are functions of compact support in Second, they are definitively not con- 
tinuous, since they are differentiable with respect to W, although they are cadlag, 
leaving some margin for development. Moreover, a simple computation shows that 
the convergence of JSf„G to JSf G holds in a weak sense, enough to estabUsh the Jf^ 
estimates needed to prove tightness, but not enough to ensure some sort of conti- 
nuity for ^„G. The key fact here is that for any test function G G the derivative 
dG/dW is a smooth function, since we have to take its usual derivative to compute 
^G. In order to benefit from this continuity, we rewrite the decomposition (I I) as 

^^«(G) = %"iG) + f4j2L n{g{ri:{x))-g{n:ix+l)))V'iGds + ^nG), 

where we have defined 

^nc- G{{x+\)/n)-G{x/n) 
W{{x+l)/n) - W{x/n)' 

In other words, V"G is a discrete approximation of dG/dW. At this point, it 
is not clear at all that rewriting the integral term in this form is of some help, 
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because now we have a big factor n in front oig{r\^{x))—g{r\^{x+l)), and this last 
difference is not small, though it should be small after time integration. Denote by 
(•, •)p the inner product in ^^(Vp). We start recalling Kipnis-Varadhan inequality, 
valid for any reversible Markov process. 

Proposition 11. For any function h: ^^Rin .if^(Vp) with J hdVp = 0, 



|2 



where the Sobolev-type norm || • „ w defined by 

||/j|p_i_„ = sup{2(/,/z)-(/, -L„/)p}, 
/ 

Ln is the generator of the process r]" and the supremum is over functions f in 

Of course, in this proposition nothing prevents us from having „ = +00, 

but estimating | |_i „ is usually part of the work when applying this inequaUty. 

Theorem 6. For any function H -.W—^R, 

2i 



( [ -TJ2 L<sKi^))-Sm^+ l)))H{x/n)dsy 

< 5t(p{p)-^ H{x/n f{W{{x+ l)/n) - W{x/n)) . 



Proof. Let us call h{ri") the term inside the integral. By Proposition 11, the expec- 
tation is bounded by 20f | \^i„. Looking at the variational formula for | ^, it 
will be good to obtain a more explicit formula for (/, — L„/)p . After some changes 
of variables and some algebra, it is not difficult to see that 

{f,-Lnf)p = I / ginix)) [/(r]^^+i)-/(r])]^p(Jr?). 

xeZ \ n I \n) 

Use now the weighted Cauchy-Schwartz inequality ah < a^P/2 + b^/2p with 
a = \/8{ri{x)) and b = V§(r](x))[/(T7^'^+i) - f{ri )] to get 

I {g{n{x+l))-g{r]{x)))f{r))dVp <^J g{r]{x))[f{n^^+')-f{n)]\{dn) 

+ ^lg{v{x))dvp. 

Notice that the last integral is bounded by (p(p)j3x/2. The factor we are esti- 
mating appears multipUed by n^/^H{x/n) in {f^,h)p. Therefore, choosing ji^ = 
\H{x/n) \{W{{x+l)/n)- W{x/n)) /2nV2 we obtain the bound 



^ xez 



xel. 
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which proves the theorem. □ 

Remember that our goal is to replace in equation (11) the term (g(T]"(^)) — 
p)^nGix/n) by (p'{p){T]^{x) - p)^G{x/n). Let us take now H = dG/dW. For 
X G M and e > 0, 

fx+e 

G{x + e)-G{x)= H{y)W{dy) 

Jx 

= H{x)(y/{x + e)-W{x)) + {H{y)-H{x))W{dy). 

Jx 



In particular, 

G(x+e)-G(x) 



W{x+e)-W{x) 



-H{x) 



1 /'-^■+^ 

< 



r \H{y)-H{x)\W{dy) 

Jx 



< swp{eH'{y);x<y <x + e}. 
From this estimate, we can perform a first replacement in (1 1). 
Lemma 11. 



UmE" 



( / -f72 i:«(gWW)-5(^;(^+i)))(v^G-//(x/«))j. 

^■^ ^ xez 



0. 



Proof. Take M > and restrict the sum above to < Mn). Then the limit holds 
by Theorem 6 and the fact that H' is right-continuous. In the set { |x| > Mn}, bound 
|V"G — H{x/n)\ by |V"G| + |//(x/«)| and use Theorem 6 again to prove that the 
expectation goes to as M — >^ oo, uniformly van. □ 

In particular, we have the decomposition 

^"(G) = %{G) + r ^2 I n{g{n^{x))-g{n^{x+ l)))H{x/n)ds + ^P{G), 
■'O " xez 

plus a rest that vanishes in ^^{F") asn^°o. We have done all this work only to 
justify the exchange of V"G by H. Now we can perform the integration by parts to 
write the integral as 

T7I I - 9{pMH{{x+ \)/n)-H{x/n))ds. 

Now we have made appear a finite approximation of H'{x/n) = .ifG(x/n). 
Since G is uniformly W-Holder continuous, using the same splitting into a box 
of size Mn used above, we can write the integral as 

f -m"L - cpip))^G{x/n)ds 

JO n I ^^y. 



plus a rest that vanishes in ^^(P") as « ^ oo. Now, finally, we have put the cadlag 
test function J^G in place of ^,G. For the sake of completeness, we repeat here 
the statement of the Boltzmann-Gibbs principle. We will give an outline of the 
proof in the next section. 



ANOMALOUS FLUCTUATIONS 



31 



Theorem 7. For any continuous function of compact support // : R M we have 



UmE" 

n— »o° 



{jl^ I {8{n':{x))-<f>{p)-<P\p){n^{x)-p))H{x/n)ds) 



2-1 

= 0. 



As a corollary, in view of Lemma 1 1 we have the decomposition 

^^,"iG) = %"{G) + (p'{p) f^;i^G)ds + ^PiG), 

Jo 

plus a rest that vanishes in ^^(P") as n — > oo. Take now a subsequence n' such 
that {^"}„ converges in distribution to some process We can assume, taking 
a further subsequence if needed, that {^"} also converges to some process 
Differently from the superdiffusive case, here the function ^G belongs to ^ and 
therefore the limiting processes ^ satisfy 

^(G) = %{G) + [' W,{^G)ds + ^t{G) 
Jo 

for any function G e As we did in the superdiffusive case, now our task is to 
prove that ^t{G) is a martingale of quadratic variation 2t(f>{p)S{G,G). By the 
invariance of Vp under the evolution of T],, the integral E"[(.^,(G))] converges to 
2t(p{p)S'{G,G). Therefore, it is enough to show that the variance of {^t{G)) goes 
to as n' ^ oo. For the superdiffusive case, it was enough to take the variance of 
the integrand and to show that it goes to 0. This is not true in the superdiffusive 
case, since the increments of the form W {x + e) — W {x) do not go to with e, at 
least not uniformly, due to the jumps of W. We state the needed result in the form 
of a lemma. 



Lemma 12. 

0, 



limE" 



( / I isKix) - cpip)) {KG^A^^WdsY 



where we have used the notation A"W = W{{x+ l)/n) — W{x/n). 

Proof. First notice that the sum under the integral in the previous expression is 
uniformly bounded by 2(f„(G,G). Therefore, we can restrict the sum to a finite 
box of arbitrary size M. By the energy estimate in Theorem 6, we that for any 
M>0, 



^"\(f I {g(i7;w)-g(i7;(x+i))}(v,"G)'A«w.) 

JO \^\^T[J„ 



2i 

< 



\x\<Mn 

<V ^ (V«G)^(A«W)^ 

' \x\<Mn 

Remember that V"G is close to H{x/n). This approximation is uniform in 
bounded intervals, in view of (12). The function H is bounded in bounded in- 
tervals, since it has a derivative in ^^(M) (which is equal to the function ^G). 
Therefore, the last sum is bounded by C{G)t{W {M) -W {-M)Y /n, which goes 
to as n — >^ oo. Fix a number / > 0. For notational convenience, assume that / is a 
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divisor of n. Required changes if it is not the case are straightforward. Repeating 
the computations above a finite number of times, we see that we can introduce a 
space average in the sum of the lemma. In other words, it is enough to show that 



goes to as n oo. But now the variance with respect to Vp of the sum inside the 
integral is equal to 



which now is of order I /I. Therefore, sending / — > oo after n and before M, we 



The rest of the proof follows like in Section 5.2, so we omit it. We just point out 
that this omitted part of the proof is simplified by the fact that we have constructed 
^ in such a way that is continuous in 

6.3. The Boltzmann-Gibbs principle. In this section we give an outline of the 
proof of Theorem 7. As we mentioned before, the proof follows closely the ideas 
introduced by Chang [7], and we adopt here the proof in Chapter 11.1 of [22]. 

As we said before, the idea is to decompose the macroscopic fluctuations of 
girisi^)) i'^to two components: one given by the projection over the conserved 
quantities and another one orthogonal to the space of conserved quantities in a 
proper sense. The intuition is simple. Consider the subspace of ^^(Vp) con- 
sisting of local functions h satisfying / hdVp = 0. An example of function in 
is ri{x) — p for any x G Z. Another example is g{ri {x)) — (p{p). A simple way to 
generate lots of local functions in is to take Lf for / local, where L is the gen- 
erator of the dynamics. What is remarkable, is that the space can be equipped 
with a norm in such a way that Jfp is the orthogonal direct sum of the subspace 
generated by {L f; f local} and {ri{x)—p;x&Z}. This decomposition is moreover 
continuous with respect to the energy estimate in Theorem 6. Fortunately, due to 
the so-called gradient condition [22] of the zero-range process, we do not need to 
prove this decomposition, but only to follow the intuition given by it. 

Therefore, the idea is to find some function / such that 



plus an error term small in some sense. Here we have written Txf{r]) = /(t^t]) 
and TxT\{z} = r]{x + z), that is, Xx is the standard shift by x. We will see that the 
space-time fluctuations of Lfij] ) are small, and that the error term will be small 
with respect to Vp , when averaged over small boxes of fixed size. 

We start with the fluctuations of Lf. Fix an intermediate scale /. We will drop 
the dependence on / from the notation, unless stated expUcitly. For simpUcity, 




L\i<M« 




complete the proof of the lemma. 



□ 



g{ri{x)) - (p{p) - (p'{p){ri{x) - p) = TxLfiri) 
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assume that the function H has a support contained in (0, 1). Define xi = il and 
define the generators L\ U'^ acting on local functions / : H — M by 



+ g(r?(^+l))[/(r?^+^'^) -/('?)]}, 



+ g(r7(x+l))[/(r]^+i'^)-/(j])]}. 

In other words, the operator L' is the generator of the dynamics restricted to the 
box A, = {xi,... ,Xi^i} and L''" is the generator of a zero-range process in the same 
box, but with uniform transition rates. 

Lemma 13. For any Lip schitz function f : Nq' — > M we have 



lim E" 



ft 1 «/'-! .2 

/ -TTi I H{^i)L''\fK)ds^ 



0. 



1=0 

Proof. By Proposition 11, the expectation is bounded by 20?||F„||^ j where we 
have used F„ as a shorthand for the sum inside the integral. Now remember the 
variational formula for ||F„||^j „: 

ll'^n||-i,«= sup {l{h,Fn)p- {h,-Lnh)p]. (13) 
/ieif2(vp) 

We have the following relations between the different Dirichlet forms: 

n/l 

Y^{h,-nh)p<{h,-L„h)p, 
1=1 

{h,-LJ^h)p < 

At this point, the best estimate we can have for the difference involving W is 
W{xi+i/n) -W{xi/n) < W{l)-W{0). But the extra factor l/n will prove to be 
useful. Now we bound each one of the terms in separately. By Cauchy- 

Schwartz inequality, 

{h,V-\„f)p<^(h,-L'\„f)„ 

In (13), this term is multiplied by H(xi/n)/n^^^. Therefore, choosing jS, = 
\H{xi/n)\{W{\) — W(0))/n^^'^, we can cancel the term involving h in (13) to get 
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the bound 



W{l)-W{0) 



n/l 



n\\-l,n ^ 



£|77(x,-/n)|(T,,./,-L''0T,,/)p. 



The whole point of introducing L''° is that this operator is translation invari- 
ant. As we mentioned before, the space averaging is crucial in order to obtain the 
Boltzmann-Gibbs principle. Due to this translation invariance, all the terms of the 
form {TxJ-, —L'-^XxiF) p are equal. We conclude that the expectation in the state- 
ment of the lemma is bounded by C/n for some constant depending only on / and 
H, which proves the lemma. Notice that in the homogeneous case (when the dy- 
namics is already translation-invariant) we would have obtained a bound of order 
C/n^ for this expectation. □ 

The point now is that time integration has already played its role, and now only 
the spatial properties of the invariant measure Vp are needed to continue. Define 
Vx{r\) = givi^)) - (P{P) - 9'{p){'nix) - P)- We have the following result: 

Proposition 12. 



where the infimum is over all the functions / : Ng' — ^ M. 

Observe that the statement of this proposition depends only on the nature of the 
invariant measure Vp and the operator L' *^ which does not carry any information 
about the environment. A detailed proof of this proposition can be found in [22]. 

Now the proof of Theorem 7 is essentially finished. The scheme is the follow- 
ing. First introduce a spatial average, substituting H{{xi + j)/n) by H{xi/n), for 
7 = 1,...,/. This introduces an error that vanishes as n ^ oo and then I —^oo due to 
the uniform continuity of H in the interval [0, 1]. Then substract to each block the 
term V'^Txif. This can be done in view of Lemma 13. And then put the expecta- 
tion inside the time integration using Cauchy-Schwartz inequality at the cost of a 
multiplicative constant t^. At this point we have arrived exactly to Proposition 12, 
which end the proof of Theorem 7. 

6.4. Annealed density fluctuations. In this section we discuss Corollary 1 . As 
we already said, we do not expect to have a functional central limit theorem for 
'3(", since we do not have a good set of test functions not depending on W. For 
continuous functions G : R ^ R of bounded support, the random variables ^^^(G) 
are well defined, but even for the generahzed Ornstein-Uhlenbeck process based on 
the usual Laplacian A, ^ is not a well defined process in the set of Radon measures, 
which is the dual space corresponding to continuous test functions of bounded 
support. Therefore, smoother test functions are needed to define properly. In 
[16], the Sobolev space J^s+d is used as space of test functions in order to construct 
In [22], the construction is carried out using the space ^ test space. 
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In any case, some smoothness of test functions is required. In our case, the very 
concept of smoothness changes with W. 

Let {«„}„, {bn}n be two sequences of random variables defined in the same 
probability space. Denote by {b„\a„) the (random) value of b„ conditioned to the 
value of a„. Now we make a simple observation. 

Lemma 14. Assume that there are two random variables a, b such that a„^ a a.s. 
aj n — > oo and such that {b„\a„) — > {b\a) in distribution as n^ o°. If the random 
variables {bn}n are W^-valued (with d < +°°), then the random vector {a„,bn) 
converges in distribution to {a,b). 

Proof. Let F be a bounded, continuous function of two variables. We want to prove 
XhatE[F{an,bn)] converges to E[F{a,b)] as n — > 0°. By the dominated convergence 
theorem, and since {bn\an) {b\a),v/e see that 

E[F{a,b„)]=E[E[F{a,b„)\a„]]^^E[E[F{a,b)]\a] =E[F{a,b)]- 

Now we just need to prove that £'[F(a„,&„)—F(a,fe„)] goes to as n — >^ 0°. Since 
F is continuous, it is uniformly continuous in compact sets. Therefore, it is enough 
to prove that {{an,b„)}„ is tight. Since {a„}„ is convergent, it is automatically 
tight. We are left with tightness of At this point we need to use the fact that 

bn G K''. More precisely, we need {b„}n to be defined in a cr-compact space. Of 
course, this is the case for M'^, since for example is the increasing union of the 
compact sets Ki = {xe Mf'; \x\ <l},lG N. 

For any fixed realization of {a„}„, {{bn\an)}n is tight. For any e > 0, to each 
reaUzation we associate a compact set K such that P((Z?„|a„) ^ K) <e. Notice that 
any K is contained in one of the sets Ki. Therefore, we have chosen an increasing 
sequence of sets A/ on the underlying probability space for which the compact set 
chosen above is contained in Ki. Taking / large enough, we obtain tightness for 

{bn}n. □ 

Now it is clear on which sense Corollary 1 holds. Take a finite collection of 
continuous functions Gi,...,G/ with bounded support and a finite collection of 
times ?!<•••<//. By Theorem 3 and an approximation procedure, the vector 
(^f'^(Gi),...,^,"''^(G,)) converges in distribution to (^,^(Gi), . . . ,^,'^(G/)). 
And now by Lemma 14, we conclude that the vector (^"(Gi), . . . ,^"(G/)) con- 
verges in distribution to the vector (^f"°(Gi), . . . , ^^^"(G/)). 

7. Current fluctuations 

In this section we prove Theorem 4. Therefore, from now on we take d = I. We 
will treat the superdiffusive and subdiffusive cases separately. We start explaining 
a generalization of the original idea of Rost and Vares [27], which works in both 
cases. 

Let us define G/(x) = (1 —x/l)~^l{x > 0), where y"*" denotes the positive part of 
y. The sequence {G/}/ converges uniformly in compacts to the Heaviside function 
Hq defined in Section 4.5. This convergence also holds in another "energy" sense. 
Of course G/ is not an admissible test function (neither in the superdiffusive case 
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nor in the subdiffusive case). In the other hand, '3^t{Gi) is well defined for any t > 0, 
since ^(•) is a white noise in ^ and therefore it can be continuously extended to 
^^{M). By a similar reasoning, the joint distributions {^^j (G/), . . . , ^?r, (G/)} are 
well defined. What is not well defined is ^^(G;) as a real- valued, right-continuous 
process; this is one of the reasons why we only obtain convergence of finite distri- 
butions in Theorem 4. Let us recall the decomposition 

^t{G)-%{G) = (p'{p) f ^,{^G)ds + ^t{G). 
Jo 

Using the energy estimate in Theorem 6 and also the explicit expression for 
the quadratic variation of ^t{G), we see that the left-hand side of this identity 
is continuous with respect to the seminorm £"{0,0)^^^. The idea is that {Gi}i 
converges to under this seminorm, which allows to define the random variable 
^^(//o) — '%{Ho) by continuity. It is important to operate carefully with this ran- 
dom variable, since the terms ^(/?o) and '3^q{Hq) are not well defined. In order to 
obtain convergence of the rescaled current to this random variable, we will see that 
the sequence {^^"(G/) — ^^"(G/)}„ converges to the rescaled current, uniformly in 
n. In particular the Umits n — oo and I ^ oo will be exchangeable. 

7.1. The superdiffusive case. In Theorem 6, we proved the energy estimate for 
the subdiffusive case. The proof for the superdiffusive case is exactly the same. 
For the reader's convenience, we restate here the energy estimate. 



Theorem 8. For any function H : 



E (^WW) - (p{p))^nH{x/n)dsy] < 5t(p{p)-'Sn{H,H). 



This estimate, together with Theorem 2 allows to obtain part i) of Theorem 4. 

For X G Z define /"(f) = n^^Jx{tn^"'). First notice that 

m - X{Gi) + %{Gi) = f4j2L isKi^)) - 9{p))^n{Ho - Gi)ds 

+ ^P{Ho-Gi). 

In principle, J^P{Hq — G/) is not well defined, since Hq — Gi is not in ^^(M). 
However, using the relations T7"(x) — T7q(x) = /"(f) and the explicit formula for 
the compensator of /^(?), this formula is easily justified. Using Theorem 8 and the 
explicit form of the quadratic variation {^"{Hq — Gi)) we obtain that 

E" [(yo"(0 - ^"(G/) + %{Gi)) '] < c{p)t£„{Ho - Gi,Ho - G/) 

for some constant c(p) not depending on /, n nor t. A simple computation shows 
that (§'„{Hq — Gi,Hq — Gi) < c/^^" for some constant c, depending only on a. We 
conclude that '3^"{Gi) — '3^q{Gi) converges to 7Q(f) as / — > oo uniformly in n. 

At this point, we need to justify the hmit '3^"{Gi) '3^t{Gi) as « ^ oo, which 
does not follows from the convergence of since G/ is not a test function. But 
this is elementary, since the sequence {^t"{-)}n is uniformly continuous under 
convergence in 
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Repeating these computations for the process we see that 

E[{^t{Gi - Gi+m) - %{Gi - Gi+m)f] < ctl'-" 

for a constant c depending only on a and p. In particular, the sequence {^(G/) — 
'3^q{Gi)}i is a Cauchy sequence and ^(//q) — ^(^o) is well defined as the limit 
of this sequence. Now we can exchange the limits in n and / to get the following 
result. 

Theorem 9. For the zero-range process with long jumps and a G (1,2), starting 
from the initial distribution Vp, 

lim -7o(m2") = ^,{Hq) - %{Ho), 

n^oo fi 

in the sense of finite-dimensional distributions, where the process ^^(//o) ~ ^(^o) 
is defined as above. 

Although we have proved convergence just for one-dimensional distributions, 
the same arguments give convergence of A:-dimensional distributions by consider- 
ing properly defined ^-dimensional vectors. The only point left in order to close 
the proof of part /) in Theorem 4 is to identify the process S^tiHo) — '3^o{Ho) as a 
fractional Brownian motion. Since the finite-dimensional distributions of ^ are 
all mean-zero and Gaussian, ^(i/o) — ^(^o) inherits this property and it is a 
mean-zero, Gaussian process. Given a Gaussian, mean-zero process with station- 
ary increments ^{t), its variance E[j^{t)'^\ identifies the process. But a simple 
scaling argument shows that the distributions of (^) t^l^^^jQiX) are the same. 
Therefore, £[(^(//o) - = ct^l"", and ^(//q) - ^(^o) is a fractional 

Brownian motion of Hurst exponent H = I /2a. 

7.2. The subdiffusive case. A simple computation shows that the choice of {G/}/ 
made in the previous section does not work in the subdiffusive case. The good 
choice is Gi{x) = (1 —W{x)/W{l))^l{x > 0). In fact, after some computations we 
see that 

^niHo — Gi,Ho — Gi) 



<on{Gm — Gi, Gm — Gi) 



wiiy 

\W{m)-W{l)\ 



W{m)W{l) 

In particular, those quantities go to as / ^ 0, uniformly in n, since W{1) ^0 
as / — oo. Repeating the arguments of the previous section, it is easy to obtain the 
following result, which is just a restatement of Theorem 4, part ii). 

Theorem 10. For the zero-range process with environment W and initial distribu- 
tion Vp, 

Um -J^itn^+^l^) = ^^{Ho) - ^^{Ho), 

n—K>° n 

in the sense of convergence of finite-dimensional distributions. 
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In this theorem, is the solution of (8) and '3^^ {Hq) - %^{Ho) is defined 
by continuity as the limit of ^^(G/) — 'S^q (Gi) when I The idea now is to 
use Lemma 14 to get the annealed result stated in part ///) of Theorem 4. Despite 
the fact that Corollary 1 is weaker than Theorem 2, in order to use Lemma 14 only 
convergence in the sense of finite-dimensional distributions is needed. We state the 
corresponding result for jQ(t) as a corollary. 

Corollary 2. Under the annealed law P (g) P", 

lim -/o(?«^+^/") = =r(0, 

where ,9^{-) has the same distribution of {Hq) — '3^q(Hq) with respect to the 
annealed law. 

Remember that any result in distribution with respect to P" has a counterpart 
for the laws ® P". In particular, in order to prove part ///) of Theorem 4, we are 
just left to prove that the process ^{t) defined above is a fractional Brownian 
motion of Hurst index o;/(2 + 2a). The self-similarity follows at once from the 
scale invariance of {t,x}x and /o(f)- Of course, for each fixed W, {t) is a 
Gaussian random variable. But it is not clear why the averaged law of ^^{t) 
is also Gaussian. The key point is that the "Gaussian character" of ^^(f) does 
not depend onW . In fact, by the construction in Section 3, the process '3/'^ has 
Gaussian finite distributions for any W. In a more formal setting, we can construct 
the processes (W, '3^^) in a probability space big enough such that there are random 
variables o{W,t), l^{W,t) such that the family is i.i.d. with common 

distribution equal to a normal distribution of mean zero and variance 1, and such 
that JT'^ it) = o{W, t) C (W, t) for any VK, ? . In other words, for a fixed time t, all the 
dependence of ,9^'^ {t) in W is encoded on its variance. From this decomposition 
plus the fact that the sum of two independent, Gaussian variables is also Gaussian, 
it is clear that ^{t) is Gaussian. Since any self-similar, Gaussian process is a 
fractional Brownian motion, we have finished the proof of Theorem 4, part Hi). 

8. The simple exclusion process with variable diffusion 

coefficient 

In this section we obtain a central limit theorem for a tagged particle in the 
simple exclusion process with variable diffusion coefficient as an application of 
the results on this article in the subdiffusive case. 

Consider a system G Z} of continuous-time, interacting random walks 

on the one-dimensional lattice Z, and let A. = {A,-; i G Z} be a given sequence of 
positive numbers. The dynamics of these particles is the following. The particle x; 
waits an exponential time of rate 2A,, at the end of which it chooses one of its two 
neighbors with equal probability. If there is no other particle Xj at the chosen site 
at that moment, the particle jumps to that site; if the site is occupied the particle 
stays where it is. In any case a new exponential time of rate 2A,- starts afresh. This 
happens independently for each particle. Notice that particles only interact through 
the so-called exclusion rule. This dynamics corresponds to a Markov process x(?) 
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defined on the state space = {x G Z / xj if / / j} and generated by tlie 
operator 



where {e,;/ G Z} is the canonical basis in Z^. The number A; is interpreted as the 
diffusion coefficient of particle Xi. We call the process x(f) the simple exclusion 
process with variable diffusion coefficient. Notice that when the sequence {A,},- 
is constant, x{t) is simply the usual simple exclusion process, but labeled: the so- 
called stirring process. Notice that the relative ordering of particles is preserved 
by the dynamics. Then, without loss of generality, we assume that < (f) 
for any / G Z and any t > 0. We call particle xo the tagged particle. We want to 
study the asymptotic behavior of xo{t). In particular, we want to obtain a central 
limit theorem for this particle. Again without loss of generality, we can assume 
that xo(0) = 0, that is, that the tagged particle is at the origin for t = 0. 

The relation between this model and the model studied in this article comes from 
the following observation. Define, for each / G Z and each t > 0, rit{i) = (t) — 
Xi(t) — l. A simple computation shows that in fact the process rit = {rjt (/) ; / G Z} is 
a zero-range process with interaction rate g{n) = > 0) and transition probability 
p{i,i — 1) = p{i — 1,0 = A,-. Remarkably, the position xo(?) of the tagged particle 
is equal to the current 7o(0 through the bond (—1,0) for rj,. Now take {A,}; as a 
realization of a sequence of i.i.d. random variables with common distribution A, 
such that A~^ is in the domain of attraction of an a -stable law, < a < 1. Then 
the process rjt is just the zero-range process with random environment defined in 
Section 4.4. 

Consider the environment process X*^(?) defined hy x^{t) =Xi{t) —XQ{t). In other 
words, ^{t) corresponds to the relative position of the particles with respect to the 
tagged particle. Knowing the relation between X°(/) and r],, it is not difficult to 
identify the invariant measures for ^{t). In fact, since the invariant measures of 
r\t are products of geometric distributions (for our particular choice of g{-)), we 
conclude that the distance between particles is a geometric distribution. Let us 
denote by fXp the invariant measure associated to x(f) corresponding to Vp, that it, 
/ip{x,+i — Xi = k} = Vp{r\{i) = A: — 1}. It is not difficult to see that the density of 
particles according to jXp is equal to 1/(1 + p). 

Just translating the results of Theorem 4 into the language of the simple exclu- 
sion process with variable diffusion coefficient, we obtain the following result. 

Theorem 11. Let x{t) be a simple exclusion process with i.i.d., a-stable diffusion 
coefficients. Assume that X^(?) is distributed according to the invariant measure 
jXp. Then, the position xoit) of the tagged particle satisfies the following central 
limit theorem: 



r^7(x) ^ A,-{ 1(X + Ci G Q.ex) [/(X + Ci) - /(x)] 




lim -;co(m2+2/« 



) = ^(0 
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in the sense of convergence of finite-dimensional distributions, when averaged over 
the diffusion coefficients. Here 3^ it) is a fractional Brownian motion of Hurst 
exponent a/(2 + 2a). 
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